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We present a ^''C-NMR study of the magnetic field driven transition to complete polarization of 
the S=l/2 antiferromagnetic Heisenberg chain system copper pyrazine dinitrate Cu(C4H4N2)(N03)2 
(CuPzN). The static local magnetization as well as the low-frequency spin dynamics, probed via 
the nuclear spin-lattice relaxation rate Tj~^, were explored from the low to the high field limit and 
at temperatures from the quantum regime (fcsT <C J) up to the classical regime (fcsT ^ J). 
The experimental data show very good agreement with quantum Monte Carlo calculations over 
the complete range of parameters investigated. Close to the critical field, as derived from static 
experiments, a pronounced maximum in T^^ is found which we interpret as the finite-temperature 
manifestation of a diverging density of zero-energy magnetic excitations at the field-driven quantum 
critical point. 



I. INTRODUCTION 

Quantum critical points (QCPs). i.e., zero- 
temperature phase transitions as a function of some 
control parameter, are likely to be at the core of 
unconventional finite-temperature behavior of many 
novel materials llJ, 131. Following the pioneering analysis 
of spin chains [3, |3| and spin ladders [Hi in external 
magnetic fields, Bose-Einstein condensation of hard-core 
bosons has been related to some phase transitions in 
quantum magnets which stem from the level crossing of 
elementary triplet excitations with the ground state at 
some critical external magnetic field B ^ Be- 

Field-induced QCPs have been under intense scrutiny 
for three and quasi-two-dimensional spin S = 1/2 dimer 
systems, i.e., TlCuClg [1, 0, B and BaCuSiaOg [3, 
for 5=1/2 ladder materials Cu2(C5Hi2N2)2Cl4 [111 
and (C.5Hi2N)2CuBr4 [13, [H, for the 5 = 1 Hal- 
dane chain Ni(C5Hi4N2)2N3(PF6) [H, for the coupled 
chain compound NiCl2-4SC(NH2)2 (DTN) 0, [3 with 
S" = 1, as well as for the effective S = 1 system 
(CH3)2CHNH3CuCl3 [13, [H. All of the latter materials 
feature a gapful zero-field state with the lowest triplet 
branch condensing as the field is increased. However, a 
similar scenario can be realized in the antiferromagnetic 
S" = 1/2 Heisenberg chain (AFHC) upon decreasing the 
field through the critical value for complete polarization 
Be- The Hamiltonian of the AFHC in an external field 
reads 

H = Y,{JS,-S,+i~gt,BB-S,), (1) 

i 

where Si are spin operators and J is the exchange energy. 
For B > Be = 2J /{guB) the lowest elementary excitation 



is a single Ising triplet which crosses the ground state at 
B ~ Be, where the system switches from complete po- 
larization into a Luttingcr liquid of dcconfined spinous. 
As for other one-dimensional (ID) systems investigated, 
i.e.; Haldanc chains and spin ladders, true gauge symme- 
try breaking for the triplet bosons will not occur at Be, 
however power-law correlations will develop, which are 
manifested in the spin-correlation functions [3, Q ■ 

Previous studies of field-driven criticality in quantum 
magnets have been focused on thermodynamic proper- 
ties. The dynamics remain a rather open issue. There- 
fore, the purpose of this letter is to shed light on the field- 
induced spin-dynamics of the AFHC. We report results of 
a nuclear magnetic resonance (NMR) study of the low fre- 
quency spin response for a wide range of parameters from 
low fields {B <C Be) to the high field limit {B > Be), 
as well as from the quantum regime (/c^T <C J) to the 
classical regime {ksT ^ J). On the one hand, the dy- 
namics are probed by the nuclear spin-lattice relaxation 
rate Tj~^ measured in the metalorganic AFHC CuPzN. 
On the other hand, the experimental data are compared 
to quantum Monte Carlo (QMC) calculations. Addition- 
ally we also investigated the static magnetic properties 
by comparison of the NMR frequency shift S with the 
magnetization calculated by QMC. We find very good 
agreement between experiment and theory in all cases. 



II. EXPERIMENTAL 

The compound CuPzN, i.e., copper pyrazine dinitrate 
Cu(C4H4N2)(N03)2, is one of the best realizations of 
the AFHC. C omp ared to oxide-based AFHC systems 
, [21I, [H, I23I I it has a small exchange-coupling con- 



2 




E 


4000 




Q. 

a. 


3000 






CO 


2000 


en 




IAIN 


1000 



90 

(5 [deg] 



FIG. 1: (a) Crystal structure of CuPzN. (b) ^^C NMR spec- 
trum for Bo = 2 T, T = 7K, and /3 = 50° in the b-c plane, (c) 
NMR center shift S of spectral lines versus /3 compared with 
a simulation of the local magnetic field [26l |. 



staiit J/ks — 10.7 K which allows experimental access 
to the parameter range of the saturation field Be ~ 14.9 
T and above 2J|. The unit cell is orthorhombic, at 
room temperature the lattice constants are a = 6.712 A, 
b = 5.142 A, and c = 11.732 A j!!], see Fig. [^a). The 



ID chains are equally spaced, their axis being parallel to 
a. The Cu(II) ions on a chain are separated by pyrazine 
rings which mediate the antiferromagnetic coupling be- 
tween the copper moments mainly via superexchange. 
Recently, three-dimensional ordering was observed at 107 



mK [27[. The ratio I Jmterchain/ Jintrachain |= 4.4x10 

indicates the highly one-dimensional character of this sys- 
tem. CuPzN has been characterized by inelastic neutron 
scattering, muon-spin relaxation, magnetothermal trans- 

port, specific heat, and magnetization measurements 

[271. [2^. I29I. [soj . All of these studies are consistent with a 
description of CuPzN in terms of the AFHC. 
Single crystals of CuPzN have been grown as described 
previously [2^. The crystal used for the measurements 
presented in this letter has the dimensions 2.6 x 4.2 x 0.8 
mm'^ and a mass of m = 9.09 mg. The ^^C nucleus in the 
pyrazine ring was used as the 7=1/2 NMR probe since 
the copper nuclei yield an experimentally very unfavor- 
able spin-spin relaxation time T2 < 10 ^s. The measure- 
ments at ^'^C were done for several fields between 2 and 
28 T and temperatures between 1.5 and 50 K. 
The low-field data were recorded in an 8 T supercon- 
ducting magnet with a modified Bruker CXP200 spec- 
trometer, applying a standard inversion-recovery spin- 
echo pulse sequence. The measurements at higher fields 
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FIG. 2: Comparison of temperature-dependent NMR shift 
5 with magnetization data calculated by QMC. For clar- 
ity, an offset has been added to the data at 13.8 T 
(+1000 ppm to NMR shift and 0.0099 (/is / T) to cal- 
culated magnetization) and 28 T (-1000 ppm and -0.0099 
( ^fl / T), respectively). The QMC errors are within symbol 
size. All solid lines are a guide to the eye. 



were done at the NHMFL, Tallahassee, in a 17 T super- 
conducting magnet and a 31 T resistive magnet with a 
home-built spectrometer, using a Carr-Purcell-Mciboom- 
Gill (CPMG) pulse-sequence for B > 20 T. 



III. NUCLEAR MAGNETIC RESONANCE 

The nuclear spin-lattice relaxation rate T^^ measures 
the spin fiuctuations at the nuclear Larmor frequency w„ 



1 



= E [Alp{q)+Al,{q)] 

q p—x,y,z 

/oo 
< Sp{q,t)Sp{~q,0) > e-'^-'dt (2) 
-00 

= ^J2^F^(l^S^{q,LJn) + F,iq)SA<},^n)]. (3) 



Here, F±{q) and Fz{q) arc the geometrical form fac- 
tors and S±{q,uj) and Szili'-^) are the dynamical struc- 
ture factors of the electronic spin system. Aa/} with 
a, (3 = x,y,z are the components of the hyperfine cou- 
pling tensor A(q). In CuPzN the ^^C nuclei are cou- 
pled to the magnetic moments of the Cu(II) electrons 
via isotropic hyperfine coupling Aisoiq), mediating only 
transverse spin fluctuations, and anisotropic dipolar cou- 
pling Aj;jp((7), mediating transverse and longitudinal spin 
fluctuations [l^]. We want to compare the experimen- 
tally and theoretically determined rates for trans- 
verse fluctuations. Therefore, the dipolar contribution 
to A(q) has to be minimized. This minimum is found 
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for the orientation j3 = 50° via a study of the angular 
dependence of the NMR shift 5 = (a;„ — jBo)/jBo , see 
Figs. [Hb) and [He) 



A. NMR frequency shift 

For a fixed orientation of the external field the NMR 
shift S is related to the magnetization M{T) via 



6{T) = A{q = Q) 



M{T) 
B ' 



(4) 



The shift 5 is compared with the calculated magnetiza- 
tion of a S=l/2 AFHC in Fig. [2 scahng the latter with 
the same factor A{q = 0) = 0.101 T / for all fields. 
A{q = 0) was determined by a least-squares fit of the 2 T 
data sets. For 2 T we find excellent agreement between 
experiment and theory, both showing a broad maximum 
around 6.5 K, reflecting the onset of antiferromagnetic 
correlations. At 13.8 T, slightly below the saturation 
field, both data sets show monotonous increase toward 
saturation with decreasing temperature. The kink near 
5 K in the experimental data is due to the proximity of 
the boiling point of liquid Helium. At 28 T, experiment 
and theory deviate below 20 K. An rf heating of the sam- 
ple can be excluded since the conditions of the CPMG 
pulse-sequence were carefully adjusted. 



B. Transverse dynamic structure factor 
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FIG. 3: Field dependence of the nuclear spin-lattice relaxation 
rate of ^^C in the critical regime. Left inset: The log-scale 
plot demonstrates the linear opening of the spin gap with 
field. Right inset: The full-scale plot highlights the maximum 
of Tj~^(_B) near the T—Q K critical field. All solid lines are a 
guide to the eye. 

used a flat default model m{Lj)^ matching the zeroth mo- 
ment of the trial spectrum. The optimal spectrum follows 
from the weighted average 



daP[a\S{q, T)\Sa±{q, i^) 



(6) 



with the probability distribution P[a\S{q^T)\ adopted 
fromRef. [s^]. 



Before turning to the T^^ data, we present our method 
of calculation for the field and temperature-dependent 
transverse dynamic structure factor S±{q,ujn)- Switch- 
ing to imaginary time r the latter reads S±{q,T) = 
\ /g°°dw5i((7,w)i4'(w,T), with a kernel i^(w,T) = e-^^+ 
Q-iP-r)uj p — 1/T. In real-space S±{q, r) can be cal- 
culated efficiently, using QMC. Following Ref. [s^l 

^g^/3"(n + l)(n-m)!m!^ 

SnP)S-ip + m))^^. , (5) 

where S±{q, r) = J2a e'''" Safi{T)/N and a, label lattice 
sites in a chain of length N. {. . ■)w refers to the Metropo- 
lis weight of an operator string of length n generated by 
the stochastic series expansion of the partition function 
(35I. [s^t. and p,m are positions in this string. 

Analytic continuation from imaginary times, i.e., 
S±{q,T), to real frequencies, i.e., S±{q,uj), is performed 
by the maximum entropy method (MaxEnt), minimizing 
the functional Q = — cea (33, El]- Here x refers 

to the covariance of the QMC data to the MaxEnt trial- 
spectrum Sa±{q,(-o). Overfitting is prevented by the en- 
tropy a = X)t^'S'a±(g,'^)ln[5a±(g,tj)/m(tj)]. We have 



C. 1/Ti: experiment versus theory 

Turning to the form factors in Eq. ([3]), Fig. [TJa) 
shows that the NMR site, i.e., the carbon nucleus, is 
located asymmetrically between two Cu(II) ions. There- 
fore on-site and next-nearest-neighbor correlations are 
included by using an effective real space form factor 
F{r) = [Fo6{r) + FiS{r — a)], where a is the lattice con- 
stant and -Fb,i parameterize the hyperfine coupling be- 
tween the nucleus and its nearest copper moments. This 
leads to a transverse relaxation rate of 

^=F^ [il + R^)S^iO,Lo) + 2RS^il,u)] l^^o, (7) 
J-i 

where S±{r = 0(1), w) are the real-space transvcrse-spin- 
correlation functions at a distance r = 0(1). 

In Fig. [3] we compare the observed NMR rate with 
the QMC results versus magnetic field in the quantum 
regime ksT < J, with T = 1.6 K. The QMC data is 
shown for i? = and a single overall scaling factor Fq 
assigned at 2 T and high temperatures. The similarity 
between experiment and theory is remarkable. For both 
we find a pronounced maximum of T-^B) at B = 13.8 
T shifting to lower fields with increasing temperature. To 
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FIG. 4: Temperature dependence of the nuclear spin-lattice 
relaxation rate of ^^C for different external fields. The solid 
QMC data lines are each polynomial fits to 50 analytic con- 
tinuations of a 128 site system and the error tube was chosen 
to contain all data points within a range of 2 a. The log-scale 
inset shows the exponential decrease in l/Ti with 1/T above 
Be 

interpret these results, we note that in the fully polarized 
state for B > Be, single magnons are exact eigenstates 
of Eq. Il]) with a dispersion of 



Ey{k) = Jcos(fc) + gfisB , 



(8) 



i?> (fc) displays a field-driven excitation gap of g^sB — J 
leading to an exponential decrease in T{'^{B) at fixed 
T and for B > B^. This can be seen for both, NMR 
and QMC, on the log-scale left inset in Fig. [31 The rates 
calculated by QMC display a broader maximum than the 
measured data, but drop with the same slope for fields 
above 16 T. We emphasize that this deviation between 
NMR data and QMC is confined to low temperatures 
T < 6 K and to a hmited range of fields 15 T < B < 17 
T which can be seen from the log-scale of the left inset. 
At B = Be the dispersion touches the zero at k = tt/2 
with a quadratic momentum dependence yielding a van- 
Hove type of critical DOS. This leads to the maximum 
in Tf \ tending to diverge as T ^ 0. For both, NMR 
experiment and QMC, the maximum in Fig. [3] occurs at 
Be ~ 13.8 T, which is slightly less than the saturation 
field of Be = 14.9 T for the magnetization. Most likely 
this downshift is a finite- temperature effect of excitations 
populating the gap. In the Luttinger liquid for B < Be 
the low-energy spinon excitations have a field-dependent 
Zmear dispersion, yielding a finite, yet reduced NMR rate. 

In Fig. m we compare Tf ^ rates observed experi- 
mentally with QMC results versus temperature for three 
fields, i.e., above, at, and below Be- As for the field 
dependence, the agreement between theory and exper- 
iment is very good. Inclusion of next-nearest neighbor 
hypcrfine couplings, i.e., R= —0.2, can slightly improve 
this agreement at high temperatures but decreases the 



agreement at low temperatures. The main result of this 
figure is the diverging NMR rate at Be which is very sug- 
gestive of critical scattering as T — > 0. As T increases, 
the van-Hove singularity in the DOS at Be is smeared 
leading to the decrease in T-f ^. For B = 28 T > Be, 
the rate has dropped by ~ 3 orders of magnitude due 
to the spin-gap and increases with temperature following 
l/Ti oc exp(— A/fcsT). The corresponding fits at 23 and 
28T, shown in the inset of Fig. 4, give A23T = 9.6r±0.6T 
and A28T = 14. 3r ± 0.9r, confirming that the gap in- 
creases with gusB — J. Finally, for B < Be we observe 
only a weak overall T dependence. In the classical regime 
fcsT J the rate is decreasing with increasing fields. 
This is indicative of an excitation spectrum dominated 
by spin-diffusion modes from q ~ 0. 



IV. CONCLUSION 

To summarize, by a complementary analysis of exper- 
iment and theory for the low-frequency spin spectrum 
of the AFHC CuPzN, as probed by the NMR Tf ^ rate 
as well as by the Knight shift, we have provided clear 
evidence for critical dynamics close to a field-induced 
QCP. Both, experiment and QMC calculations are in 
good agreement and show a pronounced maximum 
in Tj^^ in the vicinity of the saturation field, which 
tends to diverge as T — > 0. Moreover, good agreement 
between theory and experiment is also found for the 
magnetization versus temperature and field, except for a 
low-T deviation at 28 T, yet to be explored. Our findings 
may be of interest in the context of other field-induced 
QCPs as ,e.g., in TlCuClg or BaCuSisOg [l^. 
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